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Abstract 

We discuss the reconstruction of neutrino flavor ratios at astrophysical sources through the future 
neutrino-telescope measurements. Taking the ranges of neutrino mixing parameters Oij as those 
given by the current global fit, we demonstrate by a statistical method that the accuracies in 
the measurements of energy-independent ratios R = 0(1^^)/ {(j){i^e) + 0(^'r)) and S = (j){ve) / 4>{i't) 
among integrated neutrino flux should both be better than 10% in order to distinguish between the 
pion source and the muon-damped source at the 3 a level. The 10% accuracy needed for measuring 
R and S requires an improved understanding on the background atmospheric neutrino flux to 
a better than 10% level in the future. We discuss the applicability of our analysis to practical 
situations that the diffuse astrophysical neutrino flux arises from different types of sources and 
each point source has a neutrino flavor ratio varying with energies. We also discuss the effect of 
leptonic CP phase on the flavor-ratio reconstruction. 

PACS numbers: 95.85.Ry, 14.60.Pq, 95.55.Vj 
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I. INTRODUCTION 



The operation of IceCube detector [l| and the R&D effort of KM3Net [2| are important 
progresses toward a km^-sized detection capabihty in the neutrino astronomy j^. Further- 
more the radio and air-shower detectors, such as ANITA [3] and Pierre Auger detector [sl 
respectively, are also taking the data. These detectors are sensitive to neutrinos with en- 
ergies higher than those probed by IceCube and KMSNet. Finally, the radio extension of 
IceCube detector, the IceRay Q| , is also under consideration. It is expected to detect a score 
of cosmogenic neutrinos [3] per year. Motivated by the development of neutrino telescopes, 



numerous efforts were devoted to stud yin g neutrino 



neutrinos as the beam source 
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Due to the large neutrino propagation distance, the neutrino oscillation probabilities only 
depend on the mixing angles Oij and the CP phase 5 
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26l |. which make the astrophysical 



beam source favorable for extracting the above parameters, provided there are sufficient 
number of events. 

Most of the astrophysical neutrinos are believed to be produced by the decay of charged 
pion through the following chain: 7r+ yU+ + z/^ ^ e"*" + i/^ + z/e + or vr^ — /x^ + z/^ — >• 

+ + Ue + ^fi- This leads to the neutrino flux ratio (poi^e) '■ (poi^fi) '■ 4>o{^t) = 1:2:0 
at the astrophysical source where (poiiya) is the sum of and flux. Such a flux ratio 
results from an implicit assumption that the muon decays into neutrinos before it loses 
a signiflcant fraction of its energy. However, in some source the muon quickly loses its 



energy by interacting with strong magnetic flelds or with matter [27|, |28|, |29|]. Such a muon 
eventually decays into neutrinos with energies much lower than that of v^iVf^) from 7r+(7r^) 
decays. Consequently this type of source has a neutrino flavor ratio 0o(^'e) : 4>o{^fi) '■ 4>o{^t) = 
0:1:0, which is referred to as the muon-damped source. The third type of source 
emits neutrons resulting from the photo-disassociation of nuclei. As neutrons propagate to 
the Earth, Ue are produced from neutron j3 decays 30|, leading to a neutrino flavor ratio 
0o(t'e) : <Po{^fi) '■ (poi^r) = 1:0:0. Finally, neutrinos might be produced deep inside 
optically thick sources so that the flavor ratio at the source surface is signiflcantly different 
from the flavor ratio at the production point due to the oscillations 3l|]- Hence, unlike the 
previous three cases, the z/,- fraction can be signiflcant at the surface of such sources. In 
the class of sources studied by Mena et al. 311], which are referred to as the astrophysical 
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hidden sources, the neutrino flux ratio for E,j > 10^ GeV approaches to 1/3 : a : 6 with both 
a and b oscillating with the neutrino energy under the constraint a + b = 2/3. 

As mentioned before, almost all previous studies treat astrophysical neutrinos as the beam 
source for extracting neutrino mixing parameters 3^. To have a better determination of 
certain neutrino mixing parameter, for instance the atmospheric mixing angle 623 or the 
CP phase 6, a combined analysis on the terrestrially measured flavor ratios of astrophysical 
neutrinos coming from different sources, such as the pion source and the muon-damped 
source, has been considered 20|, |23|]. A natural question to ask is then how well one can 



distinguish these neutrino sources. The answer to this question depends on our knowledge of 
neutrino mixing parameters and the achievable accuracies in measuring the neutrino flavor 
ratios on the Earth such as i? = (piv^)/ {4>{i^e) + 4>{^t)) and S = (pi^Ue) / 4>{^t) ■ In this article, 
we shall provide an answer to this question with a statistical analysis. 

The possibility of measuring neutrino flavor fraction by IceCube has been discussed in 
] . It is through the measurement of muon track to shower ratio. It was demonstrated 



Ref. 



that the Ue fraction can be extracted from the above ratio by assuming flavor independence 
of the neutrino spectrum and the v^ — Vr symmetry, i.e., (t){v^i) = (pii^r)- Taking a pion source 
with E'^(j)(i'^) = 10^^ GeV cm~^s~^ 3J] and thresholds for muon and shower energies taken 



to be 100 GeV and 1 TeV respectively, the z/g fraction can be determined to an accuracy 
of 25% at IceCube for one year of data taking, or equivalently to an accuracy of 8% for a 
decade of data taking. However, the tau neutrino events are too rare to provide additional 
information on the neutrino flavor composition. The analysis in Ref. 331] as summarized 
above provides a feasibility of measuring i? in a good precision at IceCube and detectors 
with comparable capacities. In fact one may repeat the analysis in ^] and extract R and 
its associated uncertainty directly. The uncertainty in R is expected to be comparable to 
that of Uf, fraction. We note that the precisions on measuring R and S should depend on 
neutrino energies. However, for simplicity in discussions, we shall take R and S as ratios of 
integrated neutrino flux with appropriate energy thresholds for suppressing the atmospheric 
neutrino background. These ratios and the corresponding precisions, AR/R and AS/S, 
are therefore energy independent. In our analysis, we do not assume — symmetry for 
the neutrino flux measured on the Earth. We shall argue that, besides measuring R, it is 
essential to measure S in order to reconstruct the neutrino flavor ratio at the astrophysical 
source. This implies that a neutrino telescope beyond the capability of IceCube is needed 
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to study the neutrino flavor astronomy. 

It is important to understand the atmospheric neutrino background which affects the 
precisions on measuring R and S. The flux spectrum of conventional atmospheric neu- 
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361]. The measure- 



trinos which arise from pion and kaon decays is well understood 
ment on such a spectrum 37|] has reached to the energy of 10^ GeV. The prompt atmo- 
spheric neutrino flux arising from charm decays still contains large uncertainties 38, [so], 
and it has not yet been measured experimentally. The prompt atmospheric neutrino flux 
takes over the conventional one around 10^ GeV for z/g and 10^ GeV for z/^ 4l|. The 
flavor ratio of conventional atmospheric neutrino beyond TeV energies is approximately 
0c*™ (^e) '■ 0c*™ ('^/i) • 0c*™ ('^t) = 1 : 20 : 0; while the flavor ratio of prompt atmospheric 
neutrino flux is approximately 0p*™(z/e) : 0p*™(^'/x) '■ 0p*™(i^r) = 1:1: 0.1. Such flavor 
ratios differ signiflcantly from those of astrophysical neutrinos which arrive on Earth with 
0(z/^) ^ 0(^'r)- To suppress atmospheric neutrino background in the search of astrophysical 
neutrinos, energy distributions of astrophysical neutrino events and cuts on PMT hits are 
imposed [42 1. 

The paper is organized as follows. In Sec. II, we discuss properties of the probability 
matrix that links the initial neutrino flavor ratio to the ratio measured on the Earth. In Sec. 
Ill, we begin with a brief review on the current understanding of neutrino mixing angles. 
We then present the reconstructed neutrino flavor ratio at the source from the simulated 
data, which is generated by the chosen true values of the neutrino flavor ratio at the source 
and best-flt values of neutrino mixing parameters. The statistical analysis is performed with 
different measurement accuracies in R and S, as well as different ranges of neutrino mixing 
parameters. The implications of our results are discussed in Sec. IV. 



II. NEUTRINO MIXING PARAMETERS AND OSCILLATIONS OF ASTRO- 
PHYSICAL NEUTRINOS 



The neutrino flux at the astrophysical source and that detected on the Earth are related 



by 



^ 0(^e) ^ 



( P p p \ / 
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0o(i^m) 
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(1) 



where 0(1^0) is the neutrino flux measured on the Earth while (poii^a) is the neutrino flux 
at the astrophysical source, and the matrix element Pa/B is the probability of the oscillation 
1^13 ^a- The exact analytic expressions for Pq,^ are given in Eq. (lAip . It is seen that 
Pe^ = PeT- and P^^ = P^r = Ptt i^ the limit A = = i.e., ^23 = 7r/4 and = 0. In 
this case, the probability matrix P is singular with a vanishing determinant. In general, 
the determinant of this matrix remains suppressed since both A and D are expected to be 
small. For A = = D, the eigenvectors of P are given by 

1 
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with the corresponding eigenvalues 
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Aa = l, Afe = 0, Ae=-(4-3cu) 



(3) 



where = sin^2^i2. Therefore, those initial flavor ratios that differ from one another by 
a multiple of V'' shall oscillate into the same flavor ratio on the Earth. To illustrate this 
explicitly, we write the initial flux $0 at the astrophysical source as 





vV 



V2 



^'oK)-0oK))V^' + 



(4) 



where we have imposed the normalization condition (j^o^i'e) + <Po{^p.) + (poi^r) = 1- This 
normalization convention will be adopted throughout this paper. The first term on the 
right-hand side (RHS) of Eq. (jl]) can be expressed as {y/SV^ — y/QV^)/3. Hence the neutrino 
flux measured by the terrestrial neutrino telescope is 



$ = P0O = —V - —(1 - -u;)V' + ^- 



(5) 

does not 



It is seen that the vector V , with a coefficient proportional to 0o(^'^t) ~ 0o('^r;, 
appear in the terrestrially measured flux $. Hence the terrestrial measurement can not 
constrain 4>o{i'^) — (poii^r) in this case. 

The above degeneracy is lifted by either a non- vanishing ^13 {D 7^ 0) or a deviation of ^23 
from 7r/4 (A 7^ 0). To simplify our discussions, let us take D = and A 7^ 0. One can show 
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that the flux combination (1 + 4a; A/ (4 — 3a;)) 0o(i^//) — (1 — 2a; A/ (4 — 3a;)) (j)o{i^T) remains 
poorly constrained due to the suppression of det P. To demonstrate this, we observe that 

^ O A,. O/l I AN,, O/'l A\,. ^ 

„ 1 
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(6) 



-4a; 2(1 + A)cj 2(1 - A)cc; 

2(1 + A)a; (4 - a;) (1 + A2) - 2Au; (4 - u;) (1 - A^) 
\^ 2(1 - A)u; (4 - a;)(l - A^) (4 - a;)(l + A^) + 2Au; J 

ioT D — and A 0. The eigenvalues of P expanded to the second order in A are given by 
and the corresponding eigenvectors to the same order in A are 



(7) 



Y'a ^ jya 



yfb ^ js^b 



^ 2rA (1 + rA) ^ 
-l-2rA(l + rA) 
1 ^ 
\ 



-2 + 6rA 
1 - 6r A (1 - 3r A) 
1 



(8) 



with r = uj/{4 — 3uj) and A^"'^'^ the appropriate normalization factors. It is interesting to 
note that the corrections to the eigenvectors of P begin at 0{A) while the corrections to 
the corresponding eigenvalues begin at 0{A^). With the above eigenvectors, we write the 
source neutrino flux as 



$0 = 7V"1/"^ - [(1 + 4rA) 0o(t/^) - (1 - 2rA) 0o(i^r) - 2rA] N^V'^ 



+ 3 



(1 - 4rA) M^f^) + (1 - 2rA) M^r) - 3 (1 - 3rA) 



c-f/lc 



(9) 



It is easy to show that the measured flux P$o depends on V"' through the combination 



-BX'f,N^V'^ with 



B=[{1 + 4rA) 0o(z/^) - (1 - 2rA) <Po{ur) - 2rA] . 



(10) 



Clearly the flux combination (1 + 4rA) (f)o{i^ij,) — (1 — 2rA) ^ol^^r) is poorly constrained due 
to the smallness of A'j,, of the order A^. 
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TABLE I: Parameter sets chosen for our analysis 
Parameter set sin^ 6*12 sin^ 623 sin^ 6*13 6 



1 


n 09+0.02 
"-'••^^-0.02 


'-'•*^-0.06 


< 0.019 





2 


n 09+0.02 
^••J^-0.02 


'J-^'-'-o.oe 


< 0.019 





3a 


n 09+0.02 
'-'••-'^-0.02 


n 4r;+0.09 
'-'•*^-0.06 


0.016 ±0.010 





3b 


n 09+0.02 
"-'••-'^-0.02 


n 4^^+0.09 

'-'•^^-o.oe 


0.016 ±0.010 


tt/2 


3c 


n 09+0.02 

^••-'^-0.02 


n 4^^+0.09 


0.016 ±0.010 


vr 



III. STATISTICAL ANALYSIS 



To reconstruct the neutrino flavor ratio at the source with a statistical analysis, we employ 



the following best-fit values and la ranges of neutrino mixing angles 431] 

sin^ 012 = 0.32+°:°2^ sin^ 623 = 0.45^°;°^, sin^ < 0.019, (11) 

for the major part of our analysis. In the above parameter set, the best-fit value of 623 is 
smaller than 7r/4. There exist proposals to probe sin^ 623 by future atmospheric neutrino 



experiments 4^, |45| and long baseline neutrino experiments 46|. We therefore include in 
our analysis the hypothetical scenario that (sin^ 6*23) best fit = 0.55 with an error identical to 
the one associated with (sin^ 6'23)bestfit = 0.45. Finally we also consider a 613 range suggested 



by Ref. [47| where 

sin2ei3 = 0.016 ±0.010(1 a) (12) 

by a global analysis. 

In this work, we investigate uncertainties in the reconstruction of neutrino flavor ratios 
at the source for the pion source and the muon-damped source. Different choices of neutrino 
mixing parameters in our analysis are listed in Table [B Employing these mixing parameters, 
the true values of neutrino flavor ratios on the Earth and the corresponding values for R 
and S are presented in Table [Tll The true values of the neutrino flavor ratios on the Earth 
are denoted by $,r and $^ for the pion source and the muon-damped source respectively. 
They are calculated with Eq. ([T]) where P is evaluated with neutrino mixing parameters at 
their best-fit values. The flux ratios Rj^ and S"^ are obtained from $^ while R^ and S"^ are 
obtained from 



TABLE II: True values of neutrino flavor ratios on the Earth 



Parameter set = P^o,fi 



1 (024,0.37,0.39) 0.62 0.60 (0.35,0.33,0.32) 0.49 f.08 

2 (0.19,0.42,0.39) 0.71 0.51 (0.32,0.34,0.34) 0.52 0.94 
3a (0.27,0.35,0.38) 0.55 0.71 (0.36,0.33,0.31) 0.48 1.15 
3b (0.25,0.37,0.38) 0.59 0.64 (0.35,0.33,0.32) 0.49 1.07 
3c (0.23,0.40,0.37) 0.67 0.60 (0.33,0.34,0.33) 0.52 1.02 

Given a precision on measuring R, ARi/Ri, we estimate ASi/Si with two approaches. 
The first approach assumes that both ARi and ASi are dominated by the statistical errors. 



In this case, one has 



R^ 



ARi 



Si J y 1 + Ri \ Ri J ' ^"^ 

with i = 71, n [20|]. Using values of Ri and Si from Table [Tll we obtain AS-jr/S-^ = (1.1 — 
1.2)(Ai?^/i?^) and ASf^/S^ = (1.1 — 1.4)(A_R^/_R^). The second approach takes into account 
the specific complications for identifying tau neutrinos. Since tau lepton decays before it 
loses a significant fraction of its energy, tau neutrino is identified by the so-called double- 



bang or lollipop events 25 



,l33 



48( 1 ■ In IceCube or other detector with a comparable size, 



doub 
PeV 



e-bang events are observable only in a narrow energy range between 2 PeV and 20 



25 



4^ whUe the probabUity for observing a lolUpop event, though i„ctea.i„g with the 

neutrino energy, is still less than 10 ^ for = 1 EeV [3o!\. In view of these, we do not 
correlate ASi/ Si with ARi/Ri in the second approach. Rather we fix ASi/ Si while vary 
ARi/Ri for achieving the goal of distinguishing astrophysical neutrino sources. The results 
of both approaches will be presented. Before presenting the details of our analysis, we point 
out that the decays r —>■ UT-fiu^ and r — ^ Y^^lu^ each with a 18% branching ratio, produce 



extra muon events or secondary and z/^ [49|, l50|]. Cares are needed to separate these events 



from those of primary Vf. and or muons produced by the charged current interaction. 
The fitting to the neutrino flavor ratios at the source is facilitated through 




best fit 



(14) 



with i = 71, fi, (Tij. 



{ARi/Ri)Ri^ey,p, 0-5,;, oxp 



(A5'i/S'i)5'i,exp, = siri^ 9jk and a,2 



jk 



the Icr range for sj^. Here Ri^th and S'i_th are theoretical predicted values for Ri and Si 
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respectively while i?i,exp and are experimentally measured values. The values for Ri^exp 
and Si,exp are listed in Table [Tll which are generated from input true values of neutrino 
flavor ratios at the source and input true values of neutrino mixing parameters. In Ri^th 
and Si^th, the variables s'jj^ can vary between and 1 while cos 5 can vary between —1 and 



1. We note that similar functions 
the mixing angle 623 respectively 20 



lave been used for fitting the CP violation phase and 



231], assuming the source flavor ratio is known. In our 



analysis, we scan all possible neutrino flavor ratios at the source that give rise to a specific 
Xi value. Since we have taken -Rj^exp and S'j^exp as those generated by input true values of 
initial neutrino flavor ratios and neutrino mixing parameters, we have (x? )min=0 occurring 
at these input true values of parameters. Hence the boundaries for la and 3a ranges of 
initial neutrino flavor ratios are given by Axf = 2.3 and Axf = 11-8 respectively where 
^xl = xl - (xf)min = xl in our analysis. 



A. The reconstruction of initial neutrino flavor ratio by measuring R alone 

It is instructive to see how well one can determine the initial neutrino flavor ratio by 
measuring R alone. We perform such an analysis by neglecting the second term on the 
RHS of Eq. (IT^ . The la and 3a ranges for the reconstructed flavor ratios at the source are 
shown in Fig. [H For an input muon-damped source, it is seen that, with AR^/R^ = 10%, 
the reconstructed 3a range of the neutrino flavor ratio almost covers the entire physical 
region. For an input pion source with AR^^/R^^ = 10%, all possible initial neutrino flavor 
ratios are allowed at the 3a level. Clearly it is desirable to measure both R and S. 



B. The flavor reconstruction with measurements on both R and S 

In this subsection, we perform a statistical analysis with respect to simultaneous mea- 
surements of R and S. The accuracy for the measurement on R is ARi/Ri = 10% with 
i = TT, fi. Here we adopt the first approach for estimating ASi/Si while present the second 
approach in the next subsection. With the first approach, we have AS-^/S-,, = (11 — 12)% 
and ASfj,/Sfj, = (11 — 14)% depending on the parameter set chosen for calculations. 
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FIG. 1: The reconstructed ranges for the neutrino flavor ratios at the source with ARi/Ri = 10%. 
The left and right panels are results with the muon-damped source and the pion source as the 
input true source respectively. The numbers on each side of the triangle denote the flux percentage 
of a specific flavor of neutrino. The red point marks the muon-damped source $o,/i = (0, 1, 0) and 
the blue point marks the pion source <I>o,7r = (1/3,2/3,0). Gray and light gray areas respectively 
denote the la and 3a ranges for the reconstructed neutrino flavor ratios at the source. We choose 
parameter set 1 in Table |T] for this analysis. 

1. (sin^ 6*13 )best fit = 

We begin our analysis with the parameter set 1 and 2 where (sin^ ^13) best fit = and 
(sin^ 6'23)bestfit = 0.45 and 0.55 respectively. Figs. [2] and [3] show the reconstructed neutrino 
flavor ratios for an input muon-damped source and an input pion source respectively. The 
reconstructed initial flavor ratios are seen to include the region with significant Ur fractions. 
It has been shown in Sec. |ll]that the flux combination (1 -|- 4rA) 0o('^m) 2rA) (poi'^r) 

is poorly constrained due to the smallness of eigenvalue associated with V'^ (see Eq. 
and ([8])). This then leads to an extension in the reconstructed range of the initial neutrino 
flavor ratio along the V'^ direction. In the limit A = cos2^^23 = 0, i.e., sin^6'23 = 0.5, V'^ 
reduces to V'' (see Eq. (|2])) which is exactly parallel to the z/g-less side of the flavor-ratio 
triangle. The direction of V'' deviates slightly from that of V'^ in opposite ways depending 
on the sign of A. This is seen by comparing the left and right panels of both Fig. [2] and 
Fig. [31 Due to uncertainties of neutrino mixing parameters, we note that the boundaries 
for Icr and 3cr regions are not straight lines. For an input muon-damped source, the pion 
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source can be ruled out at the 3a level as shown in Fig. [2l However, the converse is not true 
as seen from Fig. [31 Finally, as shown in the right panel of Fig. |2l an astrophysical hidden 
source with $o,ah = (1/3, a, 2/3 — a) |3l|, l5l| can be ruled out at the 3a level for an input 



muon-damped source with (sin^ 6'23)bcstfit = 0.55. 



\5O0o (V,) 




FIG. 2: The reconstructed ranges for the neutrino flavor ratios for an input muon-damped source 
with AR^/Rfj, = 10% and AS^/S^ related to the former by the Poisson statistics. Gray and light 
gray areas in the left (right) panel denote the reconstructed la and 3a ranges with the parameter 
set 1 (2). 



50 00 (V,) 




4>o (Ve)'/- 



FIG. 3: The reconstructed ranges for the neutrino flavor ratios for an input pion source with 
ARjr/Rn = 10% and ASt^/S-,^ related to the former by the Poisson statistics. Gray and light gray 
areas in the left (right) panel denote the reconstructed la and 3a ranges with the parameter set 1 
(2). 



11 



2. (sin^ ) best fit > 



A non-zero ^13 introduces the CP phase contribution to every element of matrix P, except 
Pee- We study the effect of CP phase 5 on the reconstruction of neutrino flavor ratio at the 
source. We choose parameter sets 3a, 3b and 3c for performing the statistical analysis. The 
results are shown in the right panels of Figs. |l]and[5l For comparisons, we also perform the 
analysis with ^13 and ^^23 taken from the parameter set 1 and the input CP phase taken to 
be 0, 7r/2 and vr respectively. The results are shown in the left panels of Figs. H] and [51 

Left panels of Figs. H] and [5] indicate that the reconstructed ranges for initial neutrino 
flavor ratios are independent of the input CP phase for (sin^ 6'i3)bestfit = 0. The dependencies 
on the CP phase only appear in the right panels. For an input muon-damped source (see 
Fig. H]), the allowed la and 3cr ranges for initial neutrino flavor ratios are the smallest 
(denoted by red curves) for cos 5 = —1, i.e., 5 = vr. In this case, the pion source and the 



astrophysical hidden source mentioned earlier can both be ruled out at the 3(T level |51|. 
The allowed ranges become the largest (denoted by gray areas) for cos 5 = 1, i.e., 5 = 0. For 
an input pion source with different CP phases, the allowed 3cr ranges for the initial neutrino 
flavor ratio always cover the muon-damped source. 



C. Critical accuracies needed for distinguishing astrophysical sources. 

It is important to identify critical accuracies in measurement needed to distinguish be- 
tween the pion source and the muon-damped source. Choosing the parameter set 1 for the 
analysis, we present the results in Figs. [6] and [7] where two different approaches for determin- 
ing ASi/Si are used. In Fig. [6|, we determine ASi/Si by applying Poisson statistics. In the 
left panel of Fig. El which has the muon-damped source as the true source, the reconstructed 
3cr range for the neutrino flavor ratio just touches the pion source at AR^/R^ = 13% and 
AS^/S^ = 16%. In the right panel of this flgure, which has the pion source as the true 
source, the reconstructed 3a range for the neutrino flavor ratio just touches the muon- 
damped source at AP^/P^^ = 6% and AS-^/S-,, = 7%. In Fig. [71 we flx AS'^/S'^ = 25% for 
the left panel and flx AStt/Stt = 15% for the right panel. The result in the left panel is for 
AP^/P^ = 2%. We flnd that the pion source can be ruled out at the 3cr level if AP^/P^ is 
lowered to 1%. The result in the right panel is for AP^/P^ = 1.5%. If AP^/P^ is raised 
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00 (v„) 00 (Vy) 



FIG. 4: The reconstructed ranges for the neutrino flavor ratio at the source for an input muon- 
damped source with Ai?^/i?^ = 10% and AS^/Sfj_ related to the former by the Poisson statistics. 
The left panel is obtained with ^13 and 623 taken from the parameter set 1 and the input CP phase 
taken to be 0, 7r/2 and vr respectively. The right panel is obtained with the parameter sets 3a, 
3b and 3c. Light gray area, dashed blue and dashed red lines correspond to the 3a ranges for the 
reconstructed neutrino flavor ratio at the source for cos 5 = 1, cos 6 = and cos 6 = —1 respectively. 
Gray area, blue and red lines correspond to the la ranges for the reconstructed neutrino flavor 
ratio at the source for cos 6 = 1, cos 5 = and cos (5 = — 1 respectively. The effect from the CP 
phase 6 only appears in the right panel. 




00 (Vy) 00 (Vy) 



FIG. 5: The reconstructed la and 3a ranges for the neutrino flavor ratio at the source for an input 
pion source with ARj^/Rt^ = 10% and ASt^/St^ related to the former by the Poisson statistics. The 
choices of parameter sets are identical to those of Fig. HI Once more, the effect from the CP phase 
6 only appears in the right panel. 



13 



100 90 80 70 60 50 40 30 20 10 100 90 SO 70 60 50 40 30 20 10 
00 (V„) 00 (v„) 



FIG. 6: Critical accuracies needed to distinguish between the pion source and the muon-damped 
source. In the left panel where the muon-damped source is the true source, the reconstructed 3a 
range for the neutrino flavor ratio just touches the pion source at AR^/R^ = 13%. In the right 
panel where the pion source is the true source, the reconstructed 3a range for the neutrino flavor 
ratio just touches the muon-damped source at AR^^/Rt^ = 6%. We choose parameter set 1 for this 
analysis. 

to 2%, we find that the muon-damped source can not be ruled out at the 3a level. We have 
also investigated the case AS.„/S.„ = 25%. In this case the reconstructed 3a range of the 
neutrino flavor ratio covers the entire physical region unless ARt^/Rj^ is smaller than 1%. 

IV. DISCUSSION AND CONCLUSION 

The structure of the oscillation probability matrix P (singular in the limit 623 = tt/4 and 
^13 = 0) makes it difficult to constrain a flux combination approximately like the difference 
between (poi^fj.) and (pQ^u^). This then leads to an extension in the reconstructed range for 
the initial neutrino flavor ratio along the direction of V'''. 

We have illustrated the reconstruction of the neutrino flavor ratio at the source from the 
measurements of energy-independent ratios R = (piv^)/ {4>{i^e) + (/^i^r)) and 5* = 0(z/e)/0(z/^) 
among integrated neutrino flux. The ranges of neutrino mixing parameters used in this 
analysis are summarized in Eq. ( ITTi) . By just measuring R alone from either an input pion 
source or an input muon-damped source with a precision AR/R = 10%, the reconstructed 
3a range for the initial neutrino flavor ratio is almost as large as the entire physical range for 
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FIG. 7: Left panel: the reconstructed la and 3(7 ranges for the neutrino flavor ratio at the source 
for an input muon-damped source with AS^/Sfj_ = 25% and Ai?^/i?^ = 2%. Right panel: the 
reconstructed la and 3a ranges for the neutrino flavor ratio at the source for an input pion source 
with ASj^/Stt = 15% and ARt^/Rj^ = 1.5%. We choose parameter set 1 for this analysis. 

the above ratio. By measuring both R and S from an input muon-damped source, the pion 
source can be ruled out at the 3cr level for the parameter sets 1 and 2 with AR^/Rf^ = 10% 
and AS^/ related to the former by the Poisson statistics. With a pion source as the input 
true source and the choice of parameter set 1 for our analysis, the muon-damped source can 
not be ruled out at the 3a level until AR^,/ R.„ and AS-,,/ S-,, reach to 6% and 7% respectively. 
In the case (sin^ 6*13) best fit > as suggested by Ref. [47|], the CP phase 5 is seen to affect 
the reconstructed range for the neutrino flavor ratio at the source. We have also presented 
results for AS^/ and AS^r/ 5*^ fixed at 25% and 15% respectively. To distinguish the pion 
source and the muon-damped source in this case, both AR^/R^ and Ai?^/i?^ should be of 
the order 1% or smaller. 

We have also performed a statistical analysis with the errors of ^23 and 612 both reduced 
and the limit of 61^ improved to sin^ ^13 < 0.0025 (i.e., sin^ 26'i3 < 0.01). The result of 
this analysis can be best described by the modification to the left panel of Fig. [31 With 
sin^ 6'i3 < 0.0025, it is possible to rule out the muon-damped source at the Sex level for an 
input pion source by reducing the errors of both 6'23 and 612 to 70% of their original values. 

We like to point out that our analysis has been based upon the ideal scenario that the 
true astrophysical neutrino source is either a pure pion source or a pure muon-damped 
source. In practice, the neutrino flavor ratio in a single astrophysical source can depend on 
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neutrino energies such that it behaves hke the one from a pion source at the low energy 
and gradually makes a transition to the one from a muon-damped source as the neutrino 
]. Hence the reconstruction of the source flavor ratio ought to be carried 



energy mcreases 



out separately for low and high energy portions of the data. Furthermore an analysis on 
the diffuse neutrino flux is challenging since such a flux arises from astrophysical sources 
with dae.ent neutnno Bavo. .at.os. Nevertheless, the vety high epetgy part ot tlte flux 
spectrum is possibly dominated by the cosmogenic neutrino flux [52] arising from GZK 
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5J] interactions. The cosmogenic neutrino flux is a typical example of neutrino fluxes 



due to the pion source. Therefore it is sensible to reconstruct the source flavor ratio with 
respect to the highest energy part of the diffuse neutrino spectrum, provided there are 
sufficient number of events. 

The accuracy ARi/Ri = 10% (z = fi, n) frequently used in our discussions requires (9(100) 
neutrino events for each flavor. Furthermore, the above accuracy requires an improved 
understanding on the background atmospheric neutrino flux to a level better than 10% in 
the future. Taking a neutrino flux upper bound E'^cp^Ua) = 10^'' GeV cm^^s^^ (a = e, /x, r) 
derived by Waxman and Bahcall [34], we estimate by a simple re-scaling of the result in 
Ref. [55] that it takes the IceCube detector about a decade to accumulate (9(100) z/^ events. 
We stress that the above bound is for diffuse neutrino flux. The flux from individual point 
source is smaller. Hence it takes even a longer period to accumulate the same number of 
events. The IceRay jsi detector is expected to accumulate neutrino events in a much faster 
pace. However the efficiency of flavor identification in this detector still requires further 
studies. 

In summary, we have demonstrated that it is challenging to reconstruct the neutrino 
flavor ratio at the astrophysical source, requiring a lot more than a decade of data taking in a 
neutrino telescope such as IceCube for distinguishing between the pion source and the muon- 
damped source. We stress that the large uncertainty in the flavor ratios of astrophysical 
neutrinos should be taken into account as one uses these neutrinos as a beam source to 
extract the neutrino mixing parameters. 

Note added. As we were writing up this paper, we became aware of a paper by A. Esmaili 



and Y. Farzan, arXiv: 0905 . 0259 [hep-ph] , which also discusses the initial flavor composition 
of cosmic neutrinos with an approach different from ours. 
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APPENDIX A: THE EXACT OSCILLATION PROBABILITIES OF ASTRO- 
PHYSICAL NEUTRINOS 



The exact analytic expressions for the components of P are given by 



1 



1 - ^cu ) {l-D^ + D^ 



1 



1 + A) + {A - uj){l - A)D^ + 2^/iu{l - iu){l - A'^)D cos5 



uj{l - A) + (4 - cj)(l + A)L)2 _ 2^/u{l-uj){l - A2)L) cos5 
-{1- AfD^il- 



-uj[{l + A)2 + (1 - A)^D^ - (1 - A^)D\2 + 4 cos' 6)] 



1 



- -v^(l-cu)(l- A2) [(1 + A) - (1- A)/}'] Dcos6, 



p 



1 

2' 

1 

-L 

8 
1 



1-A')(1-D' + D^ 



[(1 - A')(l + 4D^ cos' 5 + Z)^) - 2(1 + A')D'] 



+ ^v^cj(l-c<;)(l- A2)A(l + D')Dcos5, 
1 



- [(1 + A')-(1 + A)'D'(1-D')] 



1 

— -UJ 



[{1 - A)' + (1 + AfD^ - (1 - A')D'(2 + 4 cos' 6)] 
+ l-Viuil - uj){l - A2) [(1 - A) - (1 + A)D'] Dcos5, 



(Al) 



where a; = sin' 2^12, A = cos2^^23, D = sin 6*13, and 6 the CP phase. 



[1] P. Berghaus [for the IceCube Collaboration], [arXiv:0812.2883l [astro-ph]. 



[2] The KMSNeT Collaboration, http://km3net.org 



[3] For a recent review, see V. S. Berezinsky, larXiv:0901.1428l [astro-ph]. 

[4] P. Gorham et al. [ANITA collaboration], |ar2avi0812J715J [astro-ph]. 

[5] T. P. A. Collaboration. larXiv:0903.3385l [astro-ph.HE]. 

[6] P. Allison et al, lar2£iYiQ9Q4J309J [astro-ph.HE]. 

17 



[7] V. S. Berezinsky and G. T. Zatsepin, Phys. Lett. 28B, 423 (1969); R. Engel, D. Seckel, and 

T. Stanev, Phys. Rev. D 64, 093010 (2001). 
[8] J. F. Beacom, N. F. Bell, D. Hooper, S. Pakvasa and T. J. Weiler, Phys. Rev. D 69, 017303 

(2004). 

[9] S. Pakvasa, Mod. Phys. Lett. A 19, 1163 (2004) [Yad. Fiz. 67, 1179 (2004)]. 
[10] M. L. Costantini and F. Vissani, Astropart. Phys. 23, 477 (2005); F. Vissani, Astropart. Phys. 

26, 310 (2006); |astro-ph /0609575 
[11] P. Bhattacharjee and N. Gupta, arXiv:hep-ph/0501191 
[12] P. D. Serpico and M. Kachelriess, Phys. Rev. Lett. 94, 211102 (2005). 
[13] P. D. Serpico, Phys. Rev. D 73, 047301 (2006). 
[14] Z. Z. Xing and S. Zhou, Phys. Rev. D 74, 013010 (2006). 
[15] W. Winter, Phys. Rev. D 74, 033015 (2006). 
[16] Z. Z. Xing, Phys. Rev. D 74, 013009 (2006). 
[17] D. Majumdar and A. Ghosal, Phys. Rev. D 75, 113004 (2007). 
[18] W. Rodejohann, JCAP 0701, 029 (2007). 
[19] D. Meloni and T. Ohlsson, Phys. Rev. D 75, 125017 (2007). 
[20] K. Blum, Y. Nir and E. Waxman, larXiv: 0706 .2070) [hep-ph]. 
[21] G. R. Hwang and S. Kim, Phys. Rev. D 78, 093008 (2008). 
[22] S. Pakvasa, W. Rodejohann and T. J. Weiler, JHEP 0802, 005 (2008). 
[23] S. Choubey, V. Niro and W. Rodejohann, Phys. Rev. D 77, 113006 (2008). 
[24] M. Maltoni and W. Winter, JHEP 0807, 064 (2008). 
[25] J. G. Learned and S. Pakvasa, Astropart. Phys. 3, 267 (1995). 

[26] H. Athar, M. Jezabek and O. Yasuda, Phys. Rev. D 62, 103007 (2000); L. Bento, P. Keranen 

and J. Maalampi, Phys. Lett. B 476, 205 (2000). 
[27] J. P. Rachen and P. Meszaros, Phys. Rev. D 58, 123005 (1998). 
[28] T. Kashti and E. Waxman, Phys. Rev. Lett. 95, 181101 (2005). 
[29] M. Kachelriess, S. Ostapchenko and R. Tomas, Phys. Rev. D 77, 023007 (2008). 
[30] L. A. Anchordoqui, H. Goldberg, F. Halzen and T. J. Weiler, Phys. Lett. B 593, 42 (2004). 
[31] O. Mena, I. Mocioiu and S. Razzaque, Phys. Rev. D 75, 063003 (2007). 

[32] For an exception, see Ref. [l^ where the authors discuss the implications of varying the initial 
neutrino flavor ratio. 



18 



[33] J. F. Beacom, N. F. Bell, D. Hooper, S. Pakvasa and T. J. Weiler, Phys. Rev. D 68, 093005 

(2003) [Erratum-ibid. D 72, 019901 (2005)]. 

[34] E. Waxman and J. N. Bahcall, Phys. Rev. D 59, 023002 (1998). 

[35] T. K. Gaisser and M. Honda, Ann. Rev. Nucl. Part. Sci. 52, 153 (2002). 

[36] G. D. Barr, S. Robbins, T. K. Gaisser and T. Stanev, Phys. Rev. D 74, 094009 (2006). 

[37] R. Abbasi et al. [IceCube Collaboration], Phys. Rev. D 79, 102005 (2009). 

[38] C. G. S. Costa, F. Halzen and C. Salles, Phys. Rev. D 66, 113002 (2002). 

[39] E. V. Bugaev, A. Misaki, V. A. Naumov, T. S. Sinegovskaya, S. I. Sinegovsky and N. Taka- 

hashi, Phys. Rev. D 58, 054001 (1998). 
[40] A. D. Martin, M. G. Ryskin and A. M. Stasto, Acta Phys. Polon. B 34, 3273 (2003). 
[41] For comparsions between conventional and prompt atmospheric neutrino flux, see, for example, 

J. F. Beacom and J. Candia, JCAP 0411, 009 (2004); H. Athar, F. F. Lee and G. L. Lin, 

Phys. Rev. D 71, 103008 (2005). 
[42] See, for example, A. Achterberg et al. [IceCube Collaboration], Phys. Rev. D 75, 102001 

(2007) ; A. Achterberg et al. [IceCube Cohaboration], Phys. Rev. D 76, 042008 (2007) 
[Erratum-ibid. D 77, 089904 (2008)]. 

[43] M. C. Gonzalez-Garcia and M. Maltoni, Phys. Rept. 460, 1 (2008) and references therein. 
[44] M. C. Gonzalez-Garcia, M. Maltoni and A. Y. Smirnov, Phys. Rev. D 70, 093005 (2004), and 

references therein. 
[45] S. Choubey and P. Roy, Phys. Rev. D 73, 013006 (2006). 

[46] P. Huber, M. Lindner, M. Rohnec, T. Schwetz and W. Winter, Phys. Rev. D 70, 073014 

(2004) . 

[47] G. L. Fogh, E. Lisi, A. Marrone, A. Palazzo and A. M. Rotunno, Phys. Rev. Lett. 101, 141801 

(2008) . 

[48] H. Athar, G. Parente and E. Zas, Phys. Rev. D 62, 093010 (2000). 

[49] F. Halzen and D. Saltzberg, Phys. Rev. Lett. 81, 4305 (1998). 

[50] J. F. Beacom, P. Crotty and E. W. Kolb, Phys. Rev. D 66, 021302(R) (2002). 

[51] It is important to claify that we focus on the energy range E^, > 10^ GeV when we refer to 
the class of sources discussed in Ref. [31]. This is to match the assumed neutrino energy range 
for the pion source and the muon-damped source. However, Mena et al. focus on the energy 
range E^, < 10^ GeV for identifying astrophysical hidden sources from the neutrino energy 

19 



spectra. 

[52] V. S. Berezinsky and G. T. Zatsepin, Phys. Lett. 28B, 423 (1969). See also R. Engel, D. Seckel 

and T. Stanev, Phys. Rev. D 64, 093010 (2001). 
[53] K. Greisen, Phys. Rev. Lett. 16, 748 (1966). 

[54] G. T. Zatepin and V. A. Kuz'min, Sov. Phys. JETP. Lett. 4, 78 (1966). 
[55] J. Ahrens et al. [IceCube Collaboration], Astropart. Phys. 20, 507 (2004). 



20 



